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ABSTRACT 



In colour-magnitude diagrams most galaxies fall in either the "blue cloud" or the 
"red sequence" , with the red sequence extending to significantly brighter magnitudes 
than the blue cloud. The bright-end of the red sequence comprises elliptical galaxies 
with boxy isophotes and luminosity profiles with shallow central cores, while fainter 
elliptical galaxies have disky isophotes and power-law inner surface-brightness profiles. 
An analysis of published data reveals that the centres of galaxies with power-law 
central surface-brightness profiles have younger stellar populations than the centres of 
cored galaxies. 

We argue that thermal evaporation of cold gas by virial-temperature gas plays 
an important role in determining these phenomena. In less massive galaxies, thermal 
evaporation is not very efficient, so significant amounts of cold gas can reach the 
galaxy centre and fill a central core with newly formed stars, consistent with the 
young stellar ages of the cusps of ellipticals with power-law surface-brightness profiles. 
In more massive galaxies, cold gas is evaporated within a dynamical time, so during 
an accretion event star formation is inhibited, and a core in the stellar density profile 
produced by dissipationless dynamics cannot be refilled. In this picture, the different 
observed properties of active galactic nuclei in higher-mass and lower-mass ellipticals 
are also explained because in the former the central supermassive black holes invariably 
accrete hot gas, while in the latter they typically accrete cold gas. 

An important consequence of our results is that at the present time there cannot 
be blue, star-forming galaxies in the most massive galactic halos, consistent with the 
observed truncation of the blue cloud at ~ L*. 

Key words: conduction - galaxies: active - galaxies: elliptical and lenticular, cD - 
galaxies: formation - galaxies: structure 



1 INTRODUCTION 

Stars form from cold gas. Disk galaxies like the Milky Way 
or the Magellanic clouds have significant quantities of cold 
gas in their disks, and this gas has sustained star forma- 
tion through most or all of the galaxies' lifetimes. Ellip- 
tical galaxies lack cold gas, and for at least several gi- 
gayears these systems have not formed significant num- 
bers of stars. Photometry of galaxies observed in the Sloan 
Digital Sky Survey shows that in a colour-magnitude dia- 
gram galaxies predominantly lie in either a "blue cloud" or 
a "red sequence", with a smaller nu mber of galaxies in a 
"green valley" between these features (|Blanton et al.l [20031 : 
iBaldrv et all |2004| : iDriver et alj 120061 ). Blue-cloud galaxies 
are forming stars, while red-sequence galaxies are not. Both 
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populations extend to faint magnitudes, but the red se- 
quence extends to significantly brighter magnitudes than the 
blue cloud. 

Within the luminous elliptical galaxies of the red se- 
quence two sub-populations can be distinguished. At very 
bright magnitudes there are objects with slightly boxy 
isophotes and luminosity profiles that become shallow 
at small radii, while at fainter magnitudes the galaxies 
mostly have disky i sophotes and powe r -law inner surface- 
brightness profiles llLauer et alJ [l995l: Faber et ail 19971 ; 



Graham. Erwin fc Asensio Ra mos 200l iTruiillo et al.lliooi : 



Lauer et al.ll2005l ; iFerrarese et al.l feOQG'). 

In this paper we argue that all these phenomena are 
reflections of the way gas at a galaxy's virial temperature 
interacts with the cold gas required for star formation. 

During an episode of star formation, only a part of 
the gas reservoir that drives the episode is converted into 
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stars; some of the residual gas is heated by supernovae to a 
temperature ~3x 10 6 K (|Larsonlll974l ; iDekel fc Silklll986h . 
Dark halos with masses 5 Merit — 10 12 Mq have potential 
wells that are too shallow to contain supernova-heated gas, 
so it flows out into intergalactic space. More massive halos 
do confine the hot gas gravitationally, with the result that 
once the halo's mass exceeds M cr i t , the density of hot gas 
in a halo builds up. As the density rises, the central cooling 
time of the gas shortens, and the rate at which gas accretes 
onto the central black hole (BH) increases. Studies of cooling 
flows suggest that the temperature sensitivity of the accre- 
tion rate onto the BH leads to an unsteady equilibrium, in 
which jets powered by ac cretion onto the BH replenish en - 
ergy radiated by the gas (jBirzan et al.ll20o4 lBinnevl f2005). 
Hence halos with M > Af cr it are filled with gas at the virial 
temperature, and the density of this gas steadily increases. 

We argue that such trapped hot gas eliminates cold gas 
by a combination thermal conduction and ablation. In the 
more massive galaxies of the red sequence, cold gas is elimi- 
nated within a dynamical time, with the result that when an 
object that contains cold gas falls into such a galaxy, the in- 
falling gas has negligible chance of reaching the galaxy's core 
before being heated to near the virial temperature. Hence, 
once a core develops in the luminosity profile of a massive 
galaxy of the red sequence, it cannot be be filled in by a cen- 
tral burst of star formation. A corollary is that these systems 
do not contain embedded disks, so their isophotes are more 
likely to be boxy than disky. 

In less luminous galaxies of the red sequence, the time 
required for cold gas to be heated is likely to be longer than 
a dynamical time, so when one of these galaxies encounters 
a gas-rich system, there is a central burst of star forma- 
tion that fills in the core in the luminosity profile that is 
produced by dissipationless dynamics. Moreover such cen- 
tral starbursts naturally account for the fact that ellipti- 
cal galaxies with power-law inner luminosity profiles have 
younger central stellar populations (see Section 15. 1[) . 

The galaxies of the blue cloud a re expected to contain 
hot gas - in the case of the Milky Wav. lSpitzeil (| 19561 ) already 
made this inference - but such gas is typically not detected 
in X-ray observations, suggesting that its density is lower 
than in red-sequence galaxies of similar mass. Thus, evapo- 
ration of cold gas is not a significant process in blue-cloud 
galaxies. However, in the last decade it has become clear 
that star-forming disk galaxies cycle their interstellar media 
through their halos severa l times over the Hubble time (e.g. 
iFraternali fe BinnevllioOrjl ). Consequently, once a potential 
well becomes so deep that it confines supernova-heated gas, 
and gas accumulates at the virial temperature, the cycling 
of the interstellar gas through the halo makes the interstel- 
lar gas vulnerable to evaporation. Star formation ceases and 
the galaxy quickly moves from the blue cloud to the red 
sequence. 

The paper is organised as follows. Section [5] defines 
the relevant parameters and presents formulae from which 
they can be calculated. Numerical values for representative 
galaxy models are presented in Section [3] When consider- 
ing accretion events one has to add the effects of encounters 
with systems that vary widely in mass; Section [4] provides 
a treatment of this problem. Sections [S] and [5] discuss the 
implications for elliptical galaxies and galaxies of the blue 
cloud, respectively. Our conclusions are in Section [7] 



Table 1. List of symbols. The shape factor is A(t) = 1 — e (oblate 
spheroidal clouds) and A(e) = (1 — e) 2 (prolate spheroidal clouds). 

isat(0 Conduction saturated for cloud size a < a sa t 

GradM Cloud evaporates for cloud size a < a ra( j 

tev(f) Evaporation time for cloud size a < a ra d 

*dyn( ? ") Time to fall to the centre from r 

a a f (r) Minimum size of cloud that can reach centre from r 

j^rad( r ) Cloud mass for semi-major axis a ra[ j: ^irA(e)p c a^ d 

M sa t(r) Cloud mass for semi-major axis a aa t: ■^TTA(e)p c a^ t 

M B {(r) Cloud mass for semi-major axis a s f. ^wA(e)p c a^ { 

M min max r >o M sf (r) 

2 TIMESCALES FOR EVAPORATION 

We want to determine the fate of a cloud of cool (T c <C 
10 6 K) gas that falls through gas at the virial tempera- 
ture Ti sm ~ 10 6 — 10 7 K. In particular, we ask whether a 
spheroidal cloud with semi-major axis length a that is on 
a sufficiently low angular-momentum orbit can reach the 
galaxy centre and form stars there. 

The motion of a cold cloud through a hot plasma is a 
complex dynamical process, involving heat conduction, ra- 
diative cooling, ram-pressure drag and ablation through the 
Kelvin-Helmholtz instability. We crudely simplify the treat- 
ment of this process by assuming that the cloud experiences 
only two opposed physical mechanisms: evaporation by ther- 
mal conduction from the hot interstellar medium (ISM) and 
condensation by radiative cooling. 

The rate at which hot gas will ablate a cold cloud must 
depend on the speed of the cloud's motion through the am- 
bient gas, and will be slowest for a stationary cloud; in this 
case a sheath of warm gas builds up around the cloud, par- 
tially insulating it from the hot ambient medium. In this 
paper we evaluate this minimum rate of ablation. We ne- 
glect the cloud's self-gravity, thus limiting ourselves to the 
case of clouds much smaller than the host galaxy. 

Depending on the physical properties of the ISM and 
on the size (and geometry) of the cloud, t he heat flux f rom 
the ISM to the cloud is either classical (|Spitzer|[l962T ) or 
saturated (|Cowie fc McKed H977I ). For given temperature 
and density of the ISM, we define the 'saturation size' 
a sat such that the heat flux saturates for a < a sat , and a 
'critical size' a ra d at which radiative cooling balances heat 
conduction: clouds bigger than a ra d cond ense the ambient 
medium, while smaller cloud s evaporate l|McKee fc Cowiel 
Il977l : iNipoti fc Binnevl liooi . and references therein). For 
clouds of size a < a ra( j we have to compare the evaporation 
time t cv with the dynamical time idyn- Only clouds for which 
tov > tdyn will survive long enough to form stars. 

Let us consider, for simplicity, a spherical galaxy for 
which we know as functions of r, the electron temperature 
Tl sm (r) and density n Cj i sm (r) in the hot atmosphere. At each 
radius we can compute the critical size a ra d{r), and the evap- 
oration time t ev (a,r) of a cloud of size a < a ra d(?")- At a 
given radius, the minimum size of a cloud for it to survive 
evaporation and end up forming stars is the "star-formation 
size" a s f(r) such that t ev (a s f,r) — tdyn(r)- Clouds bigger 
than a s { can reach the centre and there contribute to star 
formation. As the characteristic sizes are strongly dependent 
on the cloud geometry, it is more convenient to speak in 
terms of characteristic masses. Thus, for given cloud shape 
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and average mass density, we define the masses M ra< j, M sat , 
M s i of spheroidal clouds with semi-major axes a ra d, a sa t, a s f, 
respectively. To form stars in the central regions, an infalling 
cloud must survive evaporation at all radii, so only clouds 
more massive than M m i n = max,->o M s t (r) can contribute to 
central star formation. We will refer to M m i n as the "mini- 
mum cloud mass", because in a given galaxy all clouds less 
massive than M m i n will be evaporated and absorbed by the 
ISM within a dynamical time. TableQ]lists these definitions. 

2.1 Calculation of evaporation times 

Following ICowie fc Songailal (| 1977ft . we consider gas clouds 
modelled as prolate and oblate spheroids, which can rep- 
resent a wide range of geometries, from filaments to disks, 
through spherical clouds. We consider oblate and prolate 
spheroidal coordinates (u, v) related to the cylindrical co- 
ordinates (R,z) by R = A cosh u sin v, z — A sinh « cos u 
(oblate), and R — A sinhit sinu, z = A cosh u cos v (pro- 
late), where A is a length scale. The cloud is a prolate or 
oblate spheroid with surface u = uo- The semi-major and 
semi-minor axes of the cloud surface ar^3 a — A cosh uo and 
b — Asinhuo, and the cloud ellipticity is e = 1 — b/a — 
1 — tanhuo. We assume that in the interface the electron 
temperature and density are stratified with u. At the cloud 
surface (it = Mo), the electron temperature is T c <C Ti sm , 
while T — > Ti sm for u — > oo. 

In the regime of unsaturated heat conduction (a > a sat ), 
the heat flux is given by the classical Spitzer formula 

q cl = -k(T)VT, (1) 

where the thermal conductivity is 

k(T) = fK T 5/2 , (2) 

wher e k ~ 1.84 x 10" 5 (ln A)" 1 ergs" 1 cm -1 K _7/2 jSpitzerl 
1962), / < 1 is the factor by which magnetic fields 
suppress thermal conduct i on (e.g. iBinnev fc Cowiel Il98ll; 
iBohringer fc Fabianl Il989l ; iTribblel 1 1989ft . and In A is the 
Coulomb logarithm, which is only weakly dependent on n e 
and T (in the following we assume In A = 30). In the regime 
of classical heat co nduction and negli g ible r adiative cooling 
(a sat < a < a ra d)i ICowie fc Songailal (|l977ft computed an- 
alytically the evaporation rate of prolate and oblate clouds 
of cold (T c ~ 0) gas immersed in a medium of temperature 
Ti sm : the evaporative mass loss rate is 



where A(e) = 1 — e (oblate) or A(e) = (1 — e) 2 (prolate), nu.c 
is the cloud hydrogen density, and when deriving the nu- 
merical value we have used p c — 1.3ra p 7iH,c (appropriate for 
abundances Y = 0.25, X = 0.75). Combining equations ([3} 
and @, we get the evaporation time 



t cv (a,r) 



M 



25k B pca 2 C(e) 

-i5/2 
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M 12fK fJ,m p T? 



(5) 
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Myr, 



where C(e) = A(e)B(e) cosh[arctanh(l — e)], and we have 
used fi — 0.59. The time for a cloud to fall to the centre 
from r is the dynamical tim^EI 



37T 



16Gp d yn 



(6) 
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where pdvn( r ) = 3Md vn (r)/47rr 3 is the a verage mass density 
within r (e.g. IBinnev fc T rcmainc 1987). From equation © 
and from the condition t ev (a s f,r) — td yn (r), we derive the 
star-formation size 



a 3f ~ 0.32 



/ tdyn 

C(e) 10 8 yr 



l^J [Wk) kpa(7) 



As pointed out above, equation ((7} holds only if a sat < 
«sf < irad- Postponing to Appendix |X] a detailed treatment 
of the effects of saturation and radiation in the evaporation 
of spheroidal clouds, we report here the equations relating 
a sat and a ra d to the physical propertie s of the ISM. The 
satur ation size is (|Cowie fc McKegl 19771 ; ICowie fc Songailal 
1 1977ft 



a aa t = 5.0x10 



,/ 1/2 coshM / Tj 



O"0,sat</>sat 



M ~ 25k B B(e) ' (6) 

where [i is the mean gas particle mass in units of the pro- 
ton mass m p , ks is the Boltzmann constant, B(e) = ^ — 
2arctan[tanh(«o/2)] (oblate) and B(e) = |ln[tanh(uo/2)]| 
(prolate), with uo — arctanh(l — e). In the spherical limit 
e — > (uo ^> 1), equation B yields the mass loss r ate of a 
spherical cloud with radius o (|Cowie fc McKedl 19771 ). Given 
the average mass density of clouds p c , the mass of a cloud 
of semi-major axis a is 

M = ^A(e)a 3 p c ^ (4) 



where <^ sat ~ 1 and ao, sat (e) are dimensionless parameters 
(0.05 < (5o,sat < 1 for 0.99 >e>0, see Appendix |ATJ . In the 
case of significant saturation (a < a sat ), we should solve the 
equations for satu r ated h eat conduction. As pointed out by 
ICowie fc Songailal (|l977ft . the only analytic solution is for 
spherical symmetry (see equation I A6f) , while in prolate or 
oblate geometry a two-dimensional partial differential equa- 
tion must be solved numerically. However, when a < a sa t 
the evaporation time-scale computed for unsaturated ther- 
mal conduction provides a lower limit (and the correspond- 
ing star- formation size a B f an upper limit), and we shall see 
that these limits together with the analytic solution for the 
spherical case provide sufficient constraints for our purposes. 

In spherical symmetry, the critical size at which the 
cloud is radiatively stabilised is (|McKee fc Cowielll977ft 



1.9 x 10" 



if 



1/2 



(2~ism \ / ^e.ism \ , 
ICFkJ Um^J kpC - 



(9) 



1 Note the different notatio n: here a is the semi-major axis, while 
in lCowie &: Songailal Jl977l) a is the scale-length (our A). 



2 The actual time will be somewhat longer because of the drag 
force (See Appendix iBl for a discussion). 
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Figure 1. From top to bottom: electron density and temperature 
of the ISM, dynamical mass and dynamical time as functions of 
radius for the galaxy models LM (dotted lines) and HM (solid 
lines). The arrows indicate the bounding radii r max of the two 
models. 



Out of spherical symmetry the critical size cannot be com- 
puted analytically. However, we assume a ra( j as derived 
above to be approximately correct also for oblate and pro- 
late spheroids (see Appendix I A2I for a discussion). 

Summarising, at each radius in a galaxy we have three 
characteristic cloud masses: the star-formation mass M s {, 
the saturation mass M sa t, and the mass of radiatively sta- 
bilised clouds M ra( j. For given e and nn,c these masses are 
computed from equation Q, using the values of a s f, a sa t 
and o ra( j from equations (J7J), (JSJ) and Provided that 
M sa ,t < M s f < M ra< j, the maximum of M S f over all radii 
represents the minimum mass M m i n of a cool cloud that 
will survive evaporation long enough to contribute to a cen- 
tral starburst. If at any radius M B f > M ra d, our neglect of 
cooling is invalid; then we must assume M B { = M ra d, be- 
cause clouds more massive than M ra d condense the ambient 
medium. On the other hand, if M s f < M aat , M S f is just an 
upper limit to the actual value of the star- formation mass. 
At a given radius r, the characteristic masses M B f, M sat 
and M ra d scale with the suppression factor as f 3 ^ 2 , and de- 
pend on e and nn,c- The value of nn,c can be constrained 
because we expect that the cloud is either in pressure equi- 
librium with the ambient medium or overpressured (either 
by gravity or supernova-driven expansion) . In particular, we 
consider here clouds at T c ~ 10 4 K with hydrogen density 
n H , c = max[0.86n c , ism (T ism /10 4 K),lcm- 3 ]. 



3 CLOUD MASSES IN ELLIPTICALS 

Many elliptical galaxies have diffuse soft X-ray emission 
from a hot atmosphere. In this section we investigate how 
such an atmosphere affects the fate of any cold gas that 
might fall into the galaxy, for example during a merger. 



Table 2. Parameters of the model galaxies LM (low mass) and 
HM (high mass). M ga ,y. baryonic mass. Md yn ,tot : total dynamical 
mass. To: electron temperature of the ISM. n e o: electron number 
density of the ISM at r = 1 kpc. r max : bounding radius. 

Model 



LM 
HM 



Detailed radial temperature and density profiles are 
available f or only a few, mostly e xceptionally luminous 



M ga l 


Afdyn,tot 


To 


™c,0 


'"max 


M 


M 


K 


_ q 
cm ° 


kpc 


3.0 x 10 10 


2.5 x 10 11 


2.5 x 10 6 


2.7 x 10" 3 


20 


3.0 x 10 11 


3.9 x 10 12 


1.0 x 10 7 


1.1 x 10" 1 


80 



20031: iHumphrev et all 


2006: Fukazawa et all 


2006; 


David et al. 


2006): Randall. Sarazin & Irwin 


2006|; 



O'SuUivan. Sanderson fe Ponmanl 120071 ). While 



straightforward to evaluate the parameters of Table Q] 
for specific galaxies that have observationally estimated 
profiles, it is more instructive to use the data to build 
representative model galaxies and evaluate the quantities 
of Table [T] for these models. 

Our models are spherical and within some bounding 
radius r max have perfectly fiat circular-speed curves, so the 
total dynamical mass interior to r is Mdyn(r)/Mdy n ,tot = 
r/r max , where M dya ,tot = Af dyll (r max ). We take the ISM 
to be in hydrostatic equilibrium at a constant temperature 
Ti sm (r) = To. From these assumptions it follows that the 
density profile of the gas is a power law in r with a slope 
that depends on To; we adopted 



j r \ _ | n e ,o(r/kpc) 3/2 r < r max 
1 otherwise. 



We consider two models, the parameters of which are re- 
ported in Table [2] model LM represents a relatively low- 
mass, gas-poor elliptical galaxy such as NGC 3377 or 
NGC 3245, while model HM represents a high-mass, gas-rich 
elliptical galaxy such as NGC 4472 or NGC 4649. M dy „(r), 
"■c,ism(?")j *dyn(f)i an d 7ism(f) are plotted in Fig. [T]for mod- 
els LM (dotted lines) and HM (solid lines). 

For models LM and HM we computed the characteris- 
tic masses M ra d(r), M${ (r), M sat (r) and M m i n for oblate and 
prolate cloud models with ellipticity in the range < e < 
0.99. In all cases (even for thermal conduction suppression 
factor / ~ 1) we find a s f(r) <C r, where a s f(r) is the size of 
a cloud of mass M s f(r), consistent with our assumption of 
small clouds. Figure [2] shows the characteristic masses (in 
units of / 3 ' 2 Mq) as functions of radius for models HM (top 
panels) and LM (bottom panels) for three representative 
cloud shapes: gas filaments (prolate cloud with e = 0.98; 
left), spherical gas clouds (centre), and gas disks (oblate 
cloud with e = 0.98; right). Figure [3] plots M m i n as a func- 
tion of e for prolate and oblate clouds, for model HM (region 
shaded with red vertical lines) and LM (region shaded with 
blue diagonal lines). We checked these results by evaluat- 
ing the para meters of the gas rich galaxies NGC 4472 and 
NGC 4649 (llrwin fc Sarazinl Il996t Humphrey et al l 120061 : 
iRandall. Sarazin fc Irwinll2006r). and the gas-p oor galaxies 
NGC 3377 and NGC 3245 (iD avid et al J 12009 ). finding re- 
sults similar to those of the models LM and HM, respec- 
tively. Figs. [2] and [3] yield the following conclusions: 

• The star-formation mass M B f is an increasing function 
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Filaments (e=0.98) Spheres (e = 0) Disks (e=0.98) 




log r/kpc log r/kpc log r/kpc 



Figure 2. Characteristic cloud masses for models HM (top) and LM (bottom), for gas filaments (left), spherical clouds (centre), and 
gas disks (right). In each panel, the minimum mass Af m ; n (blue) is the maximum of M s { (red) over all radii smaller than r max (green); 
the dashed red line represents an upper limit of M a {] the blue shaded region represents values of M between the lower and upper limit 
of M m ; n ; the black curves correspond to M = M rat j (long dashed), and M = M sa t (short dashed). Masses are in units of / 3 / 2 Mq, where 
/ is the factor by which thermal conduction is reduced below the Spitzer value. 



of radius, partly because the dynamical time is and partly 
because of its dependence on nu,c- This indicates that cool 
clouds are more vulnerable in the outer regions of the galaxy. 
We note that in Fig. [2] the change of slope of the curves 
representing M rat j(r), M aat (r) and M s f(r) reflects the tran- 
sition between overpressured clouds (rc.H,c = 1cm -3 ; outer 
regions) and clouds in pressure equilibrium (tih.c > 1cm -3 ; 
inner regions; see Section T2.1 p . 

• In all cases M ra d is significantly higher than M s f, so 
neglecting radiative cooling in the computation of a s f is jus- 
tified. 

• Saturation becomes important at large radii, where the 
density of the ISM becomes low. So M s f at large radii is 
just an upper limit (dashed red lines in Fig. and, as a 
consequence, only upper and lower limits on M m i n can be 
derived (shaded blue areas in Fig. [2}. An exception is the 
spherical case, for which we computed M s f also in the sat- 
urated regime, by using as evaporation time M/M sa t, with 
M sa t given in equation (|A6|I . From Fig. [2] (central column 
of panels) it is apparent that M S f(r) for spherical clouds is 
almost constant in the saturated regime. As a consequence, 
Mmi n (solid symbols in Fig. [3)l lies close to the lower limit 



one would derive from the classical evaporation rate. This 
suggests that also in the case of prolate and oblate clouds 
the actual value of M m j n should be closer to the lower than 
to the upper limit. So, we expect the actual value of M m i n 
to be basically independent of r max . 

• For given model galaxy, M m i n depends on the geome- 
try of the cloud, with higher values obtained for very non- 
spherical systems (especially filaments) and the lowest val- 
ues obtained in the case of spherical clouds (see Fig. [3}. 
Thus for fixed mass, more non-spherical systems are more 
vulnerable to evaporation. 

• For fixed cloud geometry, M m i n is (in physical units) 
typically about three orders of magnitude higher in model 
HM than in model LM (Fig. [3} . More significantly, M min 
normalised to the baryonic mass of the galaxy M ga i is a fac- 
tor of ~ 100 higher in model HM than in model LM. Thus, 
proportionally bigger clouds are eliminated by evaporation 
in model HM than in model LM. 

• The actual value of Af m i n in physical units is largely 
uncertain because of the dependence on / 3//2 . In the (un- 
likely) case of Spitzer conductivity (/ ~ 1) we get M m i n ~ 
10 7 - 10 s M Q in model HM and M min ~ 10 4 - 10 5 M Q in 
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1 10 5 ; 



10 3 =- 



Figure 3. Minimum mass as a function of e for prolate (left) and 
oblate (right) clouds in units of f 3 / 2 MQ. The minimum mass 
lies in the region shaded with red vertical lines for model HM 
and in the region shaded with blue diagonal lines for model LM. 
The solid red square and blue circle indicate the values of M m i n 
obtained for spherical clouds for model HM and LM, respectively. 



.U„, ocM, / Axx 1 -™^*, 

M min //3A/„ 



(10) 



where a ~ 1.25 is the faint-end slope of the Schechter func- 
tion, and M* ~ 10 M© is the baryonic mass of an L* 
galaxy. The integrand above varies slowly below the expo- 
nential cutoff, so we can approximate the integral by the 
between its lower limit and the smaller of the upper limit 
and unity. Then the dimensionless measure of the impor- 
tance of cold accretion is 



gal 







(l 



/3M gal 

otherwise 



< mm 



(l 







Thus accretion shuts off completely for M m i n > 
/3 min (Mg a i, M*), and for M ga i > M* its importance de- 
clines faster than M~^. 

Since we require the mass function of gas-rich objects 
rather than the luminosity function of galaxies, it is not clear 
that it is appropriate to obtain a from the galaxy luminosity 
function: at early times the spectrum of cloud masses should 
be the same as the spectrum of dark-matter halos, which has 
a steeper slope a ~ 2. The flatter spectrum of the galaxy 
luminosity function is thought to arise from the difficulty 
of forming stars in low-mass halos, not from an absence of 
low-mass gas clouds. If we take a = 2, the integral above 
should be approximated by 



da; x 



ln[min(l,a:i)/xo] %o < min(l, a;i) 
otherwise 



(12) 



model LM, but for a more plausible value of the suppres- 
sion factor, / ~ 0.01 (jNipoti fc Binnevll2004l . and references 
therein), M m i n is as low as M m i n ~ 10 4 — lO 5 A/0 in model 
HM and M min ~ 10 - 100M Q in model LM. 



4 AGGREGATE GAS INFALL 

The quantity M m i n is the minimum mass an individual cloud 
must have to survive infall to the centre. To estimate the 
mass of cold gas that becomes available for central star for- 
mation, we need to integrate over the spectrum of cloud 
masses from this minimum mass to infinity. We can only 
guess at the form of the spectrum. Clearly it reflects the 
baryonic-mass spectrum of galaxies that are available to be 
captured by the galaxy under study, which might be ap- 
proximated by the baryonic-mass function of all galaxies 
truncated at the baryonic mass of the given galaxy to take 
account of the fact that an encounter with a more massive 
galaxy leads to destruction rather than growth. In reality 
this spectrum will overemphasise big clouds, both because 
no galaxy has more than a small fraction of its mass in 
cold gas, and because tidal shredding will split a galaxy into 
many clouds during a merger. Therefore a reasonable hy- 
pothesis for the cloud spectrum is the baryonic mass spec- 
trum of galaxies first truncated at the mass of the given 
galaxy and then scaled to smaller masses by a factor /3 in 
the range (0.01,0.1). 

If we identify the baryonic-mass function with the lu- 
minosity f unctio n, and approximate the latter with the 
ISchechteJ l | 19761) form, the mass of cold gas that can reach 
the centre of a galaxy of baryonic mass M ga l is 



gal 



M, 

Mgal 





In 



PM L 

otherwise 



< min ( 1, 



Again cold accretion declines in importance faster than M~J^ 
and cuts off for M min > /Jmin (M ga i, M«). 

Figure [3] shows the relative importance for the LM and 
HM models of cold accretion from a spectrum of cloud 
masses. The full curves apply if the latter has the same low- 
end slope as the mass spectrum of cold dark matter (CDM) 
halos, while the dashed curves apply if the low-end slope of 
the cloud spectrum is the same as that of the galaxy lumi- 
nosity function. In each case, two curves are plotted, showing 
the extreme cases of only spherical clouds (thick lines) and 
only filaments^ (thin lines). When 1.5 log/ — log/3 < 2, the 
relative importance of accretion in the LM and HM models 
is nearly independent of /, /3 and geometry but increases 
with a. Thus accretion will have the biggest differential im- 
pact on LM and HM systems when the low-end slope of the 
mass spectrum of clouds is steep. If 1.5 log/ — log/3 > 2, 
either because conductivity is efficient or clouds fragment 
strongly during mergers, accretion cuts off in HM systems 
and the ratio plotted in Fig. [4] shoots upwards. Accretion 
cuts off for smaller values of / is clouds are filamentary. 

In summary, the importance of cold accretion, mea- 
sured by M acc /M ga i, depends on the assumptions on 
the cloud mass spectrum and on the highly un- 
certain parameters / and /3. If 1.5 log/ — log/3 < 2, 
(M&cc I -Mg a i )lm/ {M&cc I A-fg a i ) hm will lie the range 3 - 30. If 



3 In this case we used as M n 
prolate clouds. 



the lower limit found for e = 0.98 
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Figure 4. The ratio between the quantity M acc /Af ga i in the LM 
model galaxy and the same quantity in the HM model galaxy as 
a function of the efficiency of thermal conduction / and of the 
parameter /3 defined in Section [4] (M aC c is the mass of accreted 
cold gas that can reach the galaxy centre and M ga j is the bary- 
onic mass of the galaxy). The full curves apply if the cloud mass 
spectrum has the same low-end slope as the mass spectrum of 
CDM halos, while the dashed curves apply if the slope is that 
of the galaxy luminosity function. Thick and thin curves refer to 
spherical clouds and filaments, respectively. 



processes such as the Kelvin-Helmholtz instability are effec- 
tive in breaking clouds into small fragments, 1.5 log/ — log/? 
could take on larger values and accretion would be sup- 
pressed in HM systems. 



5 IMPLICATIONS FOR ELLIPTICALS 

Figures [3] and [4] show that in an X-ray luminous galaxy a 
cold cloud has to be ~ 10 3 times more massive to reach the 
centre than in an elliptical with weak X-ray emission, and 
aggregate gas in infall available for central star formation is 
more important in the latter than in the former. Thus we 
expect cold gas and its associated star formation to have a 
much bigger impact on the centres of X-ray weak ellipticals 
than on the centres of their X-ray luminous brethren. In this 
section we discuss the nature of this impact and relate it to 
a dichotomy in the observed properties of elliptical galaxies. 

5.1 Observed properties of power- law and core 
galaxies 

On the basis of their central surface brightness (SB) pro- 
files, luminous elliptical galaxies fall into two distinct classes: 
power- law (or Sersic) galaxie s (PL Cs) , with SB profiles 



Figure 5. The age of the central stellar population versus the 
central slope of the SB profile 7 f or a subsample of the el liptical 
and lenticular galaxies studied by McDcrmid et al. (2006). Cen- 
tral ages (with error bar sl are from|McrjeTmi^et**aTl 120061 1. while 
the values of 7 are from lLauer et al .1 (|2007h . 



well represented by the ISersicI (|l968l ) law down to the cen- 
tre, and core (or core-Sersic) galaxies (CGs), with SB pro- 
files deviating from the Sersic law in the central regions 



because of the p r esence of a flat core dLauer et al. 1995; 
[ Faber et alj 1 19971: iGraham. Erwin fc Asensio Ramos! 120031 ; 
iTruiillo et al.ll2004lLauer et al.ll2005l : lFerrarese et al.ll200oT ). 
Quantitatively, the corresponding intrinsic stellar density 
profiles have inner logarithmic slope 7 > 0.5 in PLGs, and 
7 < 0.3 in CG s, with few "intermediate" galaxies, bavin 
0.3 < 7 < 0-5 jRavindranath et all l200ll : iRest et all l200ll 
lLauer et al.l l2007h . It is useful to summarise the main ob- 
served properties of PLGs and CGs. 

(i) Global galaxy properties: 

(a) CGs are typically more luminous than PLGs, with 
a dividing luminosity Ly ~ 3 x 10 10 Lvq. In particular, all 
galaxies brighter than Ly ~ 5 x 10 10 Lvq are CGs, while 
all galaxies fainter than Ly ~ 2 x 10 10 Lve a re PLGs. At 
intermediate luminositi es both types coexist (|Faber et al.l 
Il997l : lLauer et alj|2007h . 

(b) CGs have boxy isophotes and rotate slowly, 
PLGs have disky isophotes and rota te rapidly 
(|Kormendv fc Bender! [1998 : iFaber et al.lll997T ). 

(c) Almost all the brightest cluster galaxies are CGs 
(|Laine et al.ll2003l ). 

(d) All galaxies with high X-ray emission from hot 
gas (high Lx/Lb) are CGs, while ellipt icals with lower 
Lx/Lb include both CGs a nd PLGs (|Pellegrinil [19991 , 
120051 : lEUis fc O'Sulliv"anll2006l ). 

(ii) Properties of the central regions: 

(a) All ellipticals have central colour gradients, becom- 
ing redder towards the centre. CGs have weaker central 
colour gradients than do PLGs, though the differ ence is 
small and there is large scatter (|Lauer et al.ll2005l ). 

(b) CGs are rounder than PLGs at small radii. PLGs 
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typically have cen tral stellar disks, while CGs do not 
(|Lauer et al.ll2005h - 

(c) Nuclei (typically bluer than the surrounding 
gala xy) are more frequent in PLGs than in CGs 
(|Lauer et all |2005). In the Virgo cluster no galaxies 
brighter than Lb ~ 2.5 x 10 10 Lbp) have nuclei (|C6te et al.l 
120061 ; iFerrarese et alJlioOrj ). 

(d) Lower-velocity dispersion ellipticals have signifi- 
cantly younger central stella r populations than hig her- 
velocity dispersion ellipticals l|McDermid et alll2006l . and 
references therein). This should imply a correlation be- 
tween the stellar age of the inner regions of ellipticals and 
the central slope 7 of the SB profile. Combining info rma- 



during dissipationless merging, the orbiting BHs reduce 
the central star densit y by scattering stars out of the most 
tightly bound orbits (Begelman. Blandford fc Reesl Il9& 



tion on th e central age fromlMcDermid et al.l ( 20061 ) and 
on 7 from lLauer et all (|2007h 7 we considered a sample of 
20 elliptical and lenticular galaxies, which turn out to be 
nicely segregated in a diagram plotting central age versus 
7 (Fig. O, with CGs and PLGs having median central 
ages 13.2 Gyr and 3.6 Gyr, respective!}! 4 ]. 

(iii) Activity of the central black hole: 

(a) The mean ratio -L op t,nuc/iEdd of the optical nuclear 
emission and the Eddington luminosity of the central BH 
i s higher by two orders of mag nitude in PLGs than in CGs 



is higher by two orders 01 magn 
(|Capetti fc Balmaverdell2006l) . 



(b) If Pi. 4 is the global power of both core and ex- 
tended emission at 1.4 GHz CGs cover the full range 
(18 < IorPiaIWYLt T 1 <2&), while PLGs show no signifi- 
cant radio activity (|de Ruiter et al.ll2005h . 

(c) With R = Z/5GHz,mic/-I'B,nuc the ratio of the ra- 
dio to B-band luminosities of the nucleus and Rx = 
i/5GHz,nuc/iHX,nuc the ratio of radio and 2 — lOkeV 
X-ray luminosities, CGs are radio-loud (logi? ~ 3.6; 
logRx ~ —1-3), while P LGs are radio-quiet ( logfi ~ 
1.6; \ozRx ~ -3.3) JCapetti fc Balmaverdel l2006t 
iBalmaverde fc Capettill2006l ). 

5.2 Formation mechanisms of the central stellar 
profiles 

During hierarchical galaxy formation luminous galaxies 
arise from mergers of lower-luminosity systems. Lower- 
luminosity galaxies tend to be PLGs, so the natural starting 
point for discussing the origin of the CG/PLG dichotomy is 
the morphology of PLGs. That is, we assume that at early 
enough times any elliptical galaxy is a PLG, and our task is 
to explain why some systems became CGs. In this spirit one 
associates with each CG a mass deficit Md c f equal to the 
mass that would have to be added to the core to make the 
SB profile that of a PLG: typically Md e f ~ Af D hi where M D h 
is the central BH mass (|Milosavlievic et~al. 2002; Graham! 
l2004lMerritll2006h . This observation suggests the following 
popular "core-scouring" sc enario, which has been validat ed 
by N-body simulations l|Milosavlievic fc Merrittl 1200 ll fl: 



Ebisuzaki, Makino & Okumara 1991; Makino & Ebisuzak 



0; 
ki 



19961; iFaber et all Il997l; iQuinlan fc Hernquistl Il997 



Volonteri. Madau fc Haardd 120031 ). Thus in the 



scouring scenario the central slope 7 is determined by the 
merging history of the galaxy, and the mass deficit Mdef 
is naturally of the order of the masses of the participating 
BHs. The core-scouring scenario is successful in reproducing 
some of the observations, but it cannot be a complete ex- 
planatiorjf] because no purely stellar-dynamical mechanism 
can introduce a characteristic mass scale and thus explain 
why PLGs are never luminous, while CGs often are. 
Moreover, core scouring alone does not explain the origin 
of the correlation betwe en the presence of cores and strong 
diffus e X-ray emission l|Pellegrinil 2005; Ell is fc O 'Sullivan! 
2006), nor does it explain why flat cores have not formed 
in PLGs, given th at low-luminosity ellipticals such as M32 
also harbour BHs ijFaber et al.lll997l ). 



5.3 The role of thermal evaporation in 

determining the power-law/core dichotomy 

The discussion above suggests that a plausible working hy- 
pothesis is that all ellipticals at some stage in their evo- 
lution have central cores. The results of Sections [3] and 
13] and the relative youth of the central regions of PLGs 
(Fig. [5} suggest that in PLGs the core created by dissi- 
pationless dynamics has been filled in by star formation 
from cold gas that fell in during or after a merger. The ag- 
gregate mass of new stars formed will be proportional to 
the quantity A/ acc introduced in Section 3] The proposal 
requires that M acc > Mdef in PLGs and M acc < Mdef in 
CGs. We know that Mdef ~ Afbh, and Mb h is proportional 
to the total stellar mass o f the galaxy jMagorrian et ah! 
1 19981 ; IM arcom fc Hunt|[200l ). so it is reasonable to assume 
Mdef oc Mg a i. So Afacc/A'/dcf oc M acc /M g3l i, and our finding 
that Macc/Mgai is higher in hot-gas poor than in hot-gas rich 
ellipticals (Fig. [4} implies that the former are more likely to 
fill any central core. Quantitatively, this result will depend 
on the properties of the cloud mass spectrum. In Section 2] 
we concluded that when the power-law slope a of the mass 
function of clouds is similar to the faint-end slope of the 
galaxy luminosity function, the mass of gas that can reach 
the centre, normalised to the galaxy mass M ga i is at least 
a factor of 3 higher in lower-mass ellipticals than in higher- 
mass ellipticals. This difference rises to at least a factor of 
15 if a ~ 2, as it would be if the cloud mass spectrum were 
the same as that of CDM halos. For very massive galaxies 
for which Af m i n could be almost as high as /3Af* very little 
cold gas can reach the galaxy centre (see Section [4}. So, we 
expect fainter galaxies to accumulate cold gas that supports 
central star formation, while the more massive galaxies can 



4 The only outlier is NGC 4382, a core galaxy (7 = 0.01) with 
an estimated central ag e ~ 1.7 Gyr. We note that NGC 4382 has 
relatively low L X /L B llPellegrinij |2005|) and is characte rised by 
atypi cal central SB profile and inner colour gradient llLauer et al.l 
12005ft . 

5 Core s couring is not the only mechanism able to produce cored 
profiles: iNipoti, Londrillo fc Ciottil ll2006h showed that galaxies 



with SB profiles just like those of CGs can be obtained by a 
straightforward dissipationless collapse in the absence of a central 
BH. 



6 lLauer et all ||2005|) conclude "Overall, the evidence is consistent 
with the binary BH core formation mechanism, but it is not a 
strong endorsement" . 
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at most partially fill in the core created by dissipationless 
dynamics. The hypothesis that since its last major merger 
or accretion event every PLG has experienced an episode of 
central star formation from cold, rapidly rotating gas natu- 
rally explains all the observations listed in Section \5 . 1 1 under 
the headings (i) and (ii): 

(ia) It is possible for a faint, gas-poor elliptical to be a CG 
galaxy simply because it has not encountered a source of 
cold gas since its last major merger. The brightest ellipticals 
have dense virial-temperature atmospheres, so no luminous 
galaxy should be a PLG. 

(ib) The boxiness and slow rotation of CGs follows from the 
dominant role played by dissipationless dynamics in their 
formation. 

(ic) If a prerequisite for being a CG is the ability to con- 
fine hot gas, these galaxies will be, or have been, the central 
galaxies of groups and clusters. When two groups merge, 
their central galaxies will merge, thus ensuring that the cen- 
tral galaxy of the new group is also a CG. 

(id) We have seen that all X-ray luminous galaxies should 
be the central galaxies of groups or clusters. Such a galaxy 
can be a PLG only if its dense atmosphere built up after its 
last major merger. This is in principle possible but unlikely. 
When a small group falls into a rich cluster, it is possible that 
the central galaxy will be stripped of its dark halo before it 
is eaten by the central galaxy of the cluster. In this case we 
will observe a CG that is not X-ray luminous. There should 
be PLGs in all environments. 

(iia) The late burst of star formation that filled in the post- 
merger core of a PLG can be expected to be metal-rich and 
enhance the galaxy's metallicity gradient. 

(iib) The dynamical importance of rotation will increase as 
cold gas streams in, facilitating the formation of a central 
disk and/or a rapidly rotating kinematically decoupled core . 

(iic) If infall continues long after the merger i|Schweizerll999T ) 
and cold gas can reach the centre, a blue nuclear cluster can 
form. 

(iid) If stars cannot now form at the centres of CGs, the 
relative youth of the central stellar populations of PLGs is 
explained. 



and most of the energy released is mechanical, being asso- 
ciated with bipolar flows, which generate significant radio 
emission. The BH mass does not increase significantly, and 
star formation is in hibited. The classic ho t-mode accretor 
is Virgo A = M87 (|Di Matteo et a.1.1 l2003h . but FRI radio 
sources probably all belong to this class. 

Cold-mode accretors can be either radio-loud or radio- 
quiet, while hot-mode accretors are invariably radio- loud. 

In our picture all CGs must be hot-mode accretors, 
while PLGs can be either hot- or cold-mode accretors, de- 
pending on whether they happen to have encountered some 
cold gas recently. Items (iiia), (iiib) and (iiic) of Section f5. II 
are consistent with this expectation, and taken together sug- 
gest that most PLGs are in fact accreting cold gas. 

A possible expl anation of the radio-loudness of CGs 
is the suggestion of ICapetti fc Balmaverdd (|2006l ) that a 
galaxy's merging history determines both the presence of 
the core and the spin of the BH, and that the lat t er controls 
radio-loudness (see alsolWilson fc Colbert! Il995l; ICenll2007l; 
Sikora. Stawarz fc L asota 2003; IVolonteri. Sikora fc Lasota 
20071) . If BH spin were the only factor determining radio- 
loudness, the radio-loudness of a BH could only change 
slowly with time. Observations of micro-quasars show that 
a single source — regardless of whether it is powered by a 
BH or by a neutron star — often alternates short bursts 
of radio-loudness with longer pe riods of radio quiescence 
jNipoti. Blundell fc Binnevl [20051 '). In addition, the long- 
term time variability of central radio sources in galaxy clus- 
ters (once properly scaled for the BH mass) is co nsistent 
with that of microquasars (|Nipoti fc Binnevil2005D . These 
two pieces of evidence suggest that radio-loudness is more 
likely to be controlled by the accretion mode than by BH 
spin. The results of the present paper show that the ob- 
served radio properties of PLGs and CGs fit consistently in 
this scenario. 

Thus a considerable body of observational evidence 
points to a scenario in which the ability of cold gas to survive 
evaporation from the hot ISM is responsible for determining 
both the inner stellar profiles of elliptical galaxies and the 
radio emission from their central supermassive BHs. 



5.4 Mode of activity of the central black hole 

We now show that our picture also explains the connec- 
tions between galaxy morphology and active galactic nucleus 
(AGN) listed in Section \5 . 1 1 under the heading (iii). 

There are two very different modes of BH accre- 
tion and feedback, usually referred to as "cold mode" 
(or QSO mode) a nd "hot mode " (or radio-mode ) (e.g . 



Best et al.1 120051; iBinnevI 120051: 



Hardcastle. Evans fc Croston 



2007; 



Churazov et~ai1 120051 ; 
Siiacki et al.ll2007l ). In 



short, these are the main properties of the two modes: 

• Cold mode: the BH feeds from cold gas with L ~ Z/Edd, 
with most of the energy released going into photons (optical, 
UV, X-ray). The BH grows significantly in mass, and there is 
important attendant star formation. This mode is dominant 
at high redshift and leads to the correlation between BH 
and spheroid mass. The most extreme manifestations of cold 
mode accretion are luminous QSOs. 

• Hot mode: the BH feeds from hot gas at L <C Z/Edd, 



6 END OF THE BLUE CLOUD 

As outlined in the Introduction, a robust bimodality is 
observed in the properties of galaxies, which in colour- 
magnitude diagrams are nicely segregated into a red se- 
quence (extending to the brightest magnitudes) and a blue 
cloud (truncated at ~ L*). These features, together with 
the shape of bright-end of the galaxy luminosity function, 
are hard to reproduce in stan dard galaxy formation mod- 
els based on the pro posal by iRees fc Ostrikerl (Il977l ) and 
IWhite fc Reesl ll 19781 ) that galaxies form by cooling virial- 
temperature gas. 

Until recently it has been widely assumed that during 
virialization of a dark-matter halo gas is heated to the virial 
temperature and that galaxies form through the cooling of 
this hot gas. However, both analytic and numerical results 
suggest that during virialization only a fraction of the gas 
is heated to the virial temperature, and that this fraction is 
negligibl e in halos less massive than a critical mass M CT n ~ 
10 12 M Q l|Binnevlll977l ; lBirnboim fc DekelllioO^ ; iKeres et ail 
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120051 ; iDekel fc Birnboimll2006l ). Moreover, as iBinnevI {2004) 
first pointed out, observations of cooling flows suggest that 
heating by AGN-powered jets ensures that gas trapped at 
the virial temperature never cools. Hence this gas is not 
available for star formation, and galaxies can only go on 
forming stars and remain in the blue cloud to the extent that 
they can acquire cool gas. Hence the process studied above, 
of thermal evaporation of cold gas by hot, bears on the origin 
of the entire Hubble sequence as well as the morphological 
dichotomy of elliptical galaxies. 

While there is now a wide measure of agreement that in 
very massive galaxies star formation is "quenched" through 
lack of cold gas, the mechanism respons ible for cutting off 
the su pply of cold gas is controversial. IDekel fc Birnboiml 
argue that gravitationa l heating alone suffices (see 
also lKhochfar fc Ostrikerl2007h . While the efficiency of grav- 
itational heating undoubtedly rises with the clustering mass 
scale, several observations indicate that non-gravitational 
heating is important. Indeed studies of the temperature- 
luminosity correlation of diffuse X-ray sources has long pro- 
vide d strong evidence for powerful non-gravitational heat - 
ing ((Kaiser 199ll; IPonman. Sanderson fc Finoguenovl l2003). 
Other indications are the small fraction of baryons contained 
in galaxies, both in groups like the MW and in rich clus- 
ters, the significant fraction of the pr oducts of nucl eosyn- 
thesis that are in intracluster media <|Renzinj[l997h , and 
the direct observation of massive outflows f rom star-forming 
galaxies at both small a nd high redshift (jStrickland et al] 
12004 iPettini et al.ll2000h . 

The availability of non-gravitational heating is crucial 
because gravity alone will never drive gas out of dark mat- 
ter halos, and the fact that at faint magnitudes the galaxy 
luminosity function falls further and further below the mass 
function of dark halos has long been attributed to the abil- 
ity of non-gravitational heating to drive gas ou t of low-mass 
halos (|White fc Reeslll97Sl ; [Pekel fc Silk)lf986h . So it seems 
clear that energy released during star (and possibly black- 
hole) formation prevents virial temperature gas accumulat- 
ing in low-mass halos. When cool gas falls into these systems 
it is not heated to the virial temperature, so star formation 
can continue in them. 

It has been suggested that quasars are respon- 
sible for terminating star formatio n in their hosts 
l|Springel. Pi Matteo fc Hernquistll2005l ). We feel this pro- 
posal is implausible for two reasons: (i) The abundant ev- 
idence for (a) an association between the masses of black 
holes and their host bulges, (b) the coincidence between the 
peaks in the cosmic star-formation rate and the quasar lumi- 
nosity density, and (c) the rapidity of bulge formation, indi- 
cates that quasars are associated with rapid star formation, 
rather than quenching of star formation, (ii) Observations 
of disk galaxies imply that cold infall is a sustained process, 
while quasars flare up and then die. So while a quasar may 
clear its host of cold gas, it cannot be responsible for keeping 
the host clean after its death - which is the great majority 
of the host's life. In our view black holes have a vital role to 
play in quenching star formation, but their role is an indirect 
one: acting as thermostats for virial-temperature gas, which 
bears direct responsibility for quenching star formation. The 
connections between quasars and bulges arise because when 
cold gas is available on non-circular orbits (for example dur- 
ing a merger), bulge stars form rapidly on non-circular or- 



bits and black holes accrete rapidly in an environment that 
causes energy released by accretion to be rapidly degraded 
into low-energy photons th at produce only wea k feedback 
on the interstellar medium |Sazonov et al.ll2005l ). 

When a halo reaches the critical mass M cr i t two things 
happen almost simultaneously: (i) the fraction of infalling 
gas that stays below the virial temperature starts to dwin- 
dle, and (ii) the virial temperature rises above the temper- 
ature to which star and black-hole formation heat gas, so a 
hot atmosphere begins to accumulate - the temperature of 
this atmosphere is subsequently adjusted to the virial tem- 
perature by hot-mode accretion onto the black hole. We do 
not yet know how precise is the coincidence of the mass 
scales associated with the apparently independent processes 
(i) and (ii). However, it is does seem likely that this mass 
scale accounts for the upper limit to the luminosities of ob- 
jects in the blue cloud, and that thermal evaporation plays 
a significant role in setting this limit. Semi- analytic mod- 
els show that several features of the distribution of galax- 
ies in colour, magnitude and redshift can be accounted for 
if there is a criti cal halo mass above which star-formation 
is quenched (e.g. Croton et alj 120061 : ICattaneo et al.l 120061 : 
IDekel fc Birnboimll2006l) . 

A key object, that seems to be on the cusp of the 
transition across M cr i t is the edge-on spiral galaxy NGC 
5746. A Chandra observation shows that this galaxy is en- 
veloped by a nearly sp herical halo of X-ray emitting gas 
(Ras mussen et ai"1l2006l '). The authors suggest that star for- 
mation in the disk is sustained by cooling of this halo in the 
manner envisaged by I White fc Reesl (|l978l ). but it is hard to 
understand why cooling and star formation is not concen- 
trated at the centre of the bulge, where the halo is densest 
and the cooling time shortest, rather than further out in 
the disk, nor how such a nearly spherical cloud could give 
rise to a disk in circular motion. It seems much more likely 
that in NGC 5746 we see the first stages in the growth of 
the body of trapped virial-temperature gas that will shortly 
quench star formation, causing the galaxy to move onto the 
red sequence. 

The red sequence extends to faint magnitudes because 
galaxies hosted by halos that are much less massive than 
M cr it fall into halos of mass M > Af cr it and then quickly 
become red because (i) their own cold gas is stripped out 
of them by a combination of ablation, thermal conduction 
and supernova-driven flows, and (ii) their prospects of ac- 
creting fresh cold gas are small given that they are now 
orbiting inside a halo of mass M > M cr it and moving too 
fast to merge with other galaxies, which anyway have lit- 
tle cold gas to offer. iPostman et alj (|2005l ) showed that the 
density-morphology relation is driven by the growth with 
density i n the number of SO galaxies at the expense of spi - 
rals, and iBedregal. Aragon-Salamanca fc Merrifieldl (|2006t ) 
and iBarr et alj '(2007) have recently proven beyond reason- 
able doubt that spiral galaxies passively fade into SO galax- 
ies, presumably as a result of gas-starvation. Another indica- 
tion that many red-sequence galaxies orbit in more massive 
halos is the fact that they are preferentially found in dense 
environments, while blue-cloud galaxies typically r eside in 
lower-density environments (e.g. iBaldrv et al"]|2006l ). 

In summary, a galaxy will move from the blue cloud 
to the red sequence either because it sits in a primary halo 
and the mass of this has just surpassed Merit, or because 
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it starts to orbit sufficiently deeply in the potential of a 
primary halo with mass > Merit- From this hypothesis it 
should be possible to predict the luminosity functions of the 
blue cloud and the red sequence as functions of redshift. 
and the clustering properties of red-sequence galaxies. How- 
ever, to make these predictions one needs both a detailed 
picture of the clustering of dark halos and the distribution 
of virial-temper ature gas. T he former could be taken from 
the halo model l|Seriakll200ol ). but there is currently no ob- 
vious source, either theoretical or observational, from which 
to draw the latter: determining the distribution of hot gas 
observationally is hard because extended, low-density virial- 
temperature gas is a weak X-ray emitter and current X-ray 
detectors do not have useful velocity resolution. Building a 
model of the hot gas from cosmological simulations that in- 
clude baryons is problematic because of the importance for 
gas of feedback from stars and AGN. 



7 SUMMARY AND CONCLUSIONS 

Cool gas is a fundamental ingredient for galaxy forma- 
tion and evolution. Naturally, the study of cool gas in 
galaxies has been mainly focused on late-type galaxies, in 
which cool gas is most abundant. Early-type galaxies are 
rath er poor in cool gas, but th anks to recent surveys, in 
CO (ISage, Welch fc Youn3l2007l ). HI dMorganti et al.ll2006l ; 
lOosterloo et alj 2007 ) and optical lines ( Sarzi et alj l200fi( ). 
we now have a detailed picture of the properties of the cool 
gas also in these systems. In elliptical galaxies such cool gas 
is immersed in a sea of plasma at the galaxy's virial temper- 
ature, and a long-standing question is how these two phases 
coexist and interact. In previous works evaporation of cool 
gas by thermal conduction from gas in the hot phase has 
been invoked in order to e xplain the shutdown of star for- 
matio n in massive galaxies (|Binnev 2004; Dekcl fc Birnboiml 
2006) and the observed ionised-gas kin ematics relative t o 
that of the stars in early- type galaxies l|Sarzi et al.ll2007T ). 
but quantitative estimates of the importance of this process 
were lacking. 

This paper gives a quantitative study of the evaporation 
of cold gas clouds by thermal conduction from the hot ISM. 
We focused on two classes of luminous elliptical galaxies, the 
low-mass (LM) and high-mass (HM) systems, both of which 
probably have potential wells deep enough to trap gas heated 
by star formation, but which differ in the densities of their 
virial-temperature atmospheres. We have shown that in the 
LM systems accreted cold gas is much more likely to feed 
central star formation than in the HM systems. This fact 
naturally explains the observed dichotomy between PLCs 
and CGs in central surface brightness profiles, X-ray lumi- 
nosity, and activity of the central supermassive BH. 

Our exploration is based on a rather simplified analyt- 
ical model of the interaction between cold and hot gas, in 
which we account in detail only for thermal conduction and 
radiative cooling because we expect these to be the domi- 
nant physical processes. Ideally, the next step would be to 
explore the problem of cold-gas accretion and evaporation 
by the hot ISM using self-consistent three-dimensional hy- 
drodynamical simulations. However, due to resolution lim- 
its, it appears difficult to describe numerically such a multi- 
phase gas, even with state-of-the-art hydrodynamical simu- 



lations (e.g. iKaufmann et alJliooll iBrvanl |2007| . and refer- 
ences therein). Another useful step would be to add the re- 
sults obtained here to semi-analytic models of galaxy forma- 
tion |Bower et al.ll2006l : ICattaneo et al.ll2006l ; ICroton et al.l 
2006). Of course, a semi-analytic model is only as good as its 
input physics and such models are no substitute for ab-initio 
modelling from elementary physical principles. 

Stripping and ablation of cold gas from galaxies that are 
orbiting in the potential well of a system that has abundant 
virial-temperature gas undoubtedly plays a large role in de- 
termining the morphology of galaxies and the extension of 
the red sequence to faint magnitudes. This process is closely 
related to the one quantified in this paper. 

It now seems very likely that the central factor in divid- 
ing galaxies into the red sequence and the blue cloud is the 
existence of a critical halo mass above which gas accumulates 
at the virial temperature. What is still unclear is whether 
the supply of cold gas cuts off in galaxies associated with 
massive halos because infalling gas is then sho ck heated to 
the virial temperature (|Dekel fc Birnboimll2006l ). or because 
cold gas is evaporated and a blated by the hot gas (|Binnevl 
|2004| ; iNipoti fc Binnevl[2004T ). This paper shows that ther- 
mal evaporation is certainly a non-negligible process. 
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APPENDIX A: EVAPORATION OF 
SPHEROIDAL CLOUDS 

Al Saturation 

The classical heat flux (equation [1} is derived under the as- 
sumption that the mean free path of the electrons is small 
with respect to the temperature scale height T/\\7T\. When 
this condition is not satisfied, typically for very high electron 
temperature and/or very low electron density, equation Q 
does not hold and the heat flux is said to be saturated. As 
a discri minant of whethe r the h eat flux in a plasma is sat- 
urated, ICowie fc McKeei l| 19771) introduced the parameter 
a = Qd/^sat, where g c i = «(T)|VT| is the modulus of the 
classical heat flux, and the modulus of the saturated heat 
flux is 



Spheres 



g sa t = 5/ 



1/2 i 3 



(Al) 



where sat ~ 1 is a dimensionless parameter, p is the mass 
density of the gas, c s = U P/p — -J ksT / p,m p is the isother- 
mal sound speed, and P is the pre ssure. The quantity q aa , t 
in equation (|A1|I differs from that in lCowie fc McKeei (|l977i ) 
for the presence of f 1 ^ 2 , whic h we included to account for the 
effect s of magnetic fields (see ICowie fc McKedfl977t iBalbuj 
1 19861 . for a discussion). The critical value that separates the 
two regimes is a ~ 1: for a <C 1 the heat flux is unsatu- 
rated, while for it> 1, equation (|A1|) applies. In the case of 
the evaporation of a cloud modelled as a prolate or oblate 
spheroid we get 



f 1/2 a g(u,v,e) 



P, 



p{ u y 



where Pi sm = lim P(u) is the ISM pressure, 



25(/) S atPi S mCf jism a 



(A2) 



(A3) 
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( 10" K 



/ n c , ism \ 1 ( a \ 
V cm -3 / \ kpc / 



V kpc J 

is the saturation parameteiQ (|Cowie fc McKedll977l ). and 
(u, v) are the spheroidal coordinates introduced in Section[3] 
Here pi sm and c s ,i S m are the mass density and isothermal 
sound speed of the ISM, and we used pi sm ~ 1.15m p n Cj i sm 
and <^ S at ~ 1. The dimensionless function g(u,v,e) in the 
oblate case is 



g(u,v, e) 



[d(u ) ~ 8(u)] 1/5 



2[f( 



ik)) 



1 6/5 



cosh u\J sinh^ 



,(A4) 



U + COS 2 V 



where 9(x) = arctan[tanh(a;/2)] , while in the prolate case 



7 Ou r definition of on differs from that in l|Cowie fc Songailal 
Il977l) by a facto r cosh up, so that in th e limit e — > we recover 
the definition of lCowie &: McKe 1 lll977l) for a spherical cloud of 
radius a. 




0.01 



Figure Al. Top: do, sat ^ a function of the ellipticity e of prolate 
(left panel) and oblate (right panel) clouds. 



g(u,v,e) = 



[ 7 ( Wo )-7(«)] 1/5 



^(uo)] 6 ^ 5 sinh u\J sinh 2 u + sin 2 v 



(A5) 



ln[tanh(a;/2)], and we recall that uo 



where 7(1) 
arctanh(l — e). 

The term Pi am /P(«), appearing in equation ()A2[I . is 
expected to be of the order of unity or smaller, because 
dense clouds tend to be overpressured with respect to 
their environments, and the pressure is approximately con- 
stant throughout the fl ow if the Mach number is not large 
(Cowiefc McKeei Il977l ). While in the case of a spherical 
cloud of radius a the co ndition for saturation J max(a) > 1] 
is equivalent to ao > 1 l|Cowie fc McKeei 1 19771) , for signifi- 
cantly non-spherical prolate and oblate clouds this transition 
occurs at a critical value of ao <C 1, because the geometry 
implies stronger temperature gradients near the equatorial 
plane in the oblate case and near the symmetry axis in the 
prolate case. This is reflected in the behaviour of the func- 
tion g(u,v,e) appearing in equation ()A2|) : for e — > 1 the 
function g(u,v,e), at fixed v, has a maximum that signif- 
icantly exceeds unity, especially for v — > 7r/2 (oblate) and 
v — > (prolate). However, even for ao significantly larger 
than the critical value, the heat flux is overestimated only 
in a small region around v = tv/2 (oblate) or v = (pro- 
late), which gives a negligible contribution to the heat flow. 
If the surface of the cloud is S^ot , and the surface in which 
saturation occurs is S sa ,t, the fraction of the heat flow that 
is overestimated is C = f q c i ■ d 2 S/ f q c i ■ d 2 S, where 

d 2 S is the area element. Integration gives £ = 1 — cosw sa t 
(oblate), ( — cosw sa t (prolate), where ii sa t is the value of v 
at which the transition between classical and saturated heat 
flux occurs. As long as £ <C 1, the unsaturated computation 
still yields a good approximation of the actual heat flow. 
So, we define a reference value of the saturation parameter 
ao.sat (depending on e) such that ( < 1/3 for ao < ao, S at- 
a , S at = a , S at(e)/ao,sat(0) as a function of e is plotted in 
Fig. IA1I The mass loss rate M sa t for saturated heat flux 
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can be co mputed analytically on ly for spherically symmet- 
ric clouds l|Cowie fc McKeei[l977l) . being 

M sat = 4Trf 1/2 a 2 p isul c B>iBnl <l) Bat F(a ), (A6) 

with F(ao) — 2.73af/ 8 if sat ~ 1, where we included the 
factor f 1 / 2 to account for the effects of magnetic fields. 

A2 Radiation 

iMcKee fc Cowiel (| 19T7f > showed that a spherically symmet- 
ric cold cloud immersed in a hot (Tj 8m > 10 5 K) ambient 
medium is radiatively stabilised when oo ~ 0.027/^> sa t- For 
given physical properties of the ambient medium, this con- 
dition is satisfied when the cloud radius is a ra( j given in 
equat ion j^jl , which is approxi mate ly valid also in the oblate 
case rtCowie fc Songailalll977l ). In iNipoti fc Binnevl (|2004) 
we computed the critical length ? cf it of radiatively stabilised 
filaments of cool gas immersed in a hot medium. So, for very 
elongated prolate clouds (e — > 1) an estimate of the critical 
size is 

IcxiX „ „_ ,1/2 1/2 I Ti sm \ 7 / 4 / n c i sm \ 1 . . 

a rad ~ — c 0.07/ A ) ( — ) kpc, (A7) 

where 0.08< A(T ism ) <0.37 for 10 6 < T ism / K < fO 8 
|Nipoti fc Binnevl 12004 ). This results has been obtained in 
cylindrical symmetry for clouds with electron temperature 
T c ~ 10 4 K, assuming constant pressure in the interface. An 
approximated estimate of the evaporation time for cylindri- 
cal clouds of length / < l crit is tev.cvi ~ (l/lcrit) 2 t coo i, where 
t coo i is the cooling time (|Nipoti fc Binnevll2004[ ) . We verified 
that in the temperature range 10 6 — fO 7 K t cv , C yi ~ t cv , 
where t EV is the evaporation time of a cold prolate cloud 
with e ~ 0.98, and equations <)A7[) for filaments and (|9} for 
spheres yield to values of a ra d within a factor of two. This 
suggests that we can safely assume that equation @ gives 
an approximated estimate of a ra( j also for prolate clouds. 



t ev (M B f,r) — td yn {M B {,r). We found that the values of 
-Mmin = max r >o M s f (r) obtained considering equation (|B2|) 
are larger than those obtained in the absence of the drag 
force by 10 — 15 per cent in the model galaxies LM and HM. 
We conclude that neglecting the effect of drag implies only a 
modest underestimation of the value of the minimum mass 

This paper has been typeset from a Tp^X/ JOTjnjX file prepared 
by the author. 



APPENDIX B: DRAG FORCE 

Beside gravity, a cloud of cold gas infalling with velocity v 
through a galaxy's ISM experiences the drag force 

F drag = -^CSpismV 2 -, (Bl) 

2 v 

where S is the cloud cross-section, pi sm is the density of 
the ambient medium at the position of the cloud, and C is 
the dra g coefficient, which depen ds on the Reynolds number 
R (e.g. lLandau fc Lifshi j Il959h . In our applications R <g 
10 s , and we crudely assume C(R) = 24/ R for R < 48 and 
C(R) = 0.5 for R > 48. 

To estimate the effect of the drag force on the dynamics 
of the clouds we consider here the simple case of a spherical 
cloud of mass M infalling radially in the model galaxies LM 
and HM described in Section [3] Integrating the equation of 
motion 

r - - GM ;r (r) + ^C(r)S(r) PiBm (r)r 2 (B2) 

between r and with v(r) — 0, we obtain a dynamical 
time tdy n (M, r). We note that at fixed mass, the cross- 
section depends on r because the cloud density does (see 
Section 12.111 . The star-formation mass Af s f is such that 



